Part B

getting started

libraries

library(sf)
library(tmap)

read analysis results from Part A

rivers_south <- st_read("rivers_with_site_information.gpkg")

## Reading layer `rivers_with_site_information' from data source 
##   `/Users/mosesvelano/Desktop/OneDrive/GISC101 Semester 2/Research Project/rivers_with_site_information.gpkg' 
##   using driver `GPKG'
## Simple feature collection with 580 features and 18 fields
## Geometry type: MULTILINESTRING
## Dimension:     XY
## Bounding box:  xmin: 1090125 ymin: 4762408 xmax: 1697138 ymax: 5503850
## Projected CRS: NZGD2000 / New Zealand Transverse Mercator 2000

average_mi_south <- st_read("average_mi_south.gpkg")

## Reading layer `average_mi_south' from data source 
##   `/Users/mosesvelano/Desktop/OneDrive/GISC101 Semester 2/Research Project/average_mi_south.gpkg' 
##   using driver `GPKG'
## Simple feature collection with 243 features and 20 fields
## Geometry type: POINT
## Dimension:     XY
## Bounding box:  xmin: 1175248 ymin: 4834855 xmax: 1680898 ymax: 5446506
## Projected CRS: NZGD2000 / New Zealand Transverse Mercator 2000

ta_south_v2 <- st_read("TA_south1.gpkg")

## Reading layer `TA_south1' from data source 
##   `/Users/mosesvelano/Desktop/OneDrive/GISC101 Semester 2/Research Project/TA_south1.gpkg' 
##   using driver `GPKG'
## Simple feature collection with 22 features and 8 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: 1089970 ymin: 4747987 xmax: 1721500 ymax: 5516917
## Projected CRS: NZGD2000 / New Zealand Transverse Mercator 2000

map 1

## tmap mode set to interactive viewing

map 2

tm_shape(ta_south_v2) + tm_polygons(col = "mean_of_rivers", title = "mean of MCI by TAs")